pub fn new_birthday_probability(n: u32) -> f64 {
    let mut ans=1_f64;
    for i in 0..n{
        ans*=(365-i) as f64/365_f64;
    }
    1_f64-ans
}
